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Abstract 

This paper presents guidelines Tor developing 3 
multivariable centralized automatic flight control 
cyst-™ (AFCS) for a twin lift helicopter system 
(TLHS). Singular value ideas are used to formulate 
performance and stability robustness specifications. 
A Linear-Quadratic-Gaussian with Loop Transfer 
Recovery (LQG/LTR) design is obtained and evaluated. 


I. Introduction 

1. 1 Background 

During recent years the develcoment of heavy lift 
helicopters has leveled off. Users lave found it 
difficult to justify the high costs associated with 
such a development, particularly when there are much 
more economical alternatives. Twin lift is such an 
.alternative. It. aHovs the use of already existing 
helicopters to caoucr. missions wnicn would 
otherwise oe impossible. 

This paper illustrates how the LQG/LTR design 
methodology [5] [61, coupled with singular value ideas 
[71 [8l can be used to develop an automatic 
centralized feedback controller for 3 TLHS. 

• 2 Outline of Paper 

Because of the Inherent coupling in TLHS's, singiQ 
input-single output (SISO) design techniques can not 
easily be used to systematically develop a 
centralized AfCS. The paramount contribution of this 
paper i.: to present guidelines fer designing a ‘real 
worir multi-input multi-output (dlMO) AfCS Tor a 
TLHS. 


•This research v/as supported by the NASA Ames 
and Langley Research Centers under grant 
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The paper begins with a description of a model for 
the TLHS to be studied. A set of output variables are 
then selected. (These are the variables to bo actually 
commandad by the pilot). 

The control problem is then formulated by discussing 
the structure of the final AFCS and giving 
performance and stability robustness specifications, j 
Next, the LQG/LTR design methodology is described 
and used to obtain an AFCS. The AFCS is then 
evaluated by discussing the properties of the final 
loop, sensitivity, and closed loop trarsrer function 
matrix singular values. The characteristics of 
selected system transient responses to typical pilot 
reference commands are also discussed. 

Finally, the flexibility of the methodology to meet 
prescribed pilot specifications is discussed. 

2. Modeling a Twin Lift Helicopter System — 
2.1 Cesrrfrtico of sustem 

For simplicity this paper only focusses on controlling _ 
the longitudinal dynamics of a TLHS near hovering 
trim. The- configuration to be studied, called the 
longitudinal confi guratio n, is shown in Fig. 1. ^ 
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Fig. 1: Twin Lift Helicopter Systcn (TLHS); 
Longitudinal Configuration. 
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The longitudinal configuration consists of two 
helicopters, two tethers, a spreader bar, two load 
cables, and a payload. In studying this configuration 
it is assumed that the lateral dynamics of the 
helicopters and load-bar assembly decouple from the 
planar (longitudinal) dynamics of interest; i.e. the 
system motion Is restricted to tlw vertical 
(longitudinal) plana and hence only vertical 
translation, horizontal translation, and pitching are 
‘allowed*. 

The helicopters being modeled are Sikorsky 'JH-6QA 
Blacknawks. They are assumed to be identical. The 
lead helicopter is referred to as the master 
helicopter, and the trail helicopter as the slave. Each 
tether is assumed to have a fixed length. H. The 
helicopter-tether attachment points are assumed to 
lie a fixed distance, h, below the center of gravities 
(c.g.'s) of the helicopters. Each tether is corrected to 
one end of the spreader bar which is assumed to be 
rigid and have a fixed length. L The payload is 
suspended a fixed distance, 2. below the spreader bar 
eg. via two fixed length cables. For simplicity the 
tethers, as well as the load cabies, are assumed to 
have no compliance and to always be in tension. 

22 Variable definitions and Slate Space Model 
A seven degree of freedom linear model 111 is used to 
characterize the rigid body dynamics of the system. 
The key degrees of freedom are described as follows: 
three degrees of freedom per helicopter (vertical 
translation, horizontal translation, and pitching), and ' 
ona degree cf freedom for the lead which is modeled 
as a point mass (constrained pendular motion). 

Since this study restricts the system motion to the 
longitudinal plane, only four controls (two per 
helicopter) are relevant. These are the collective 
controls 3nd the cyclic controls. The collectives 
control the up/down motion of the helicopters. The 
cyclics contrat their fore/aft motion as well as their 
pitching motion. To fully exploit the symmetry of the 
*eoual tether configuration* being addressed, the 
following seven degrees of freedom [ll and four 
controls are defined (Fig. 2) : 

Ex n 0^ ♦ x^/Z : Average horlz. helic. coord. (I) 

Ez « (z m ♦ zJ/2 •* Average vert, helic coord. (2) 

10 a (e m ♦ 9 $ )/2 : Average helic. pitch attitude (3) 

«lc- (3 lem* ?ics^2 : Average cyclic control (4) 


E0 C a (0cm* 0(3 V2 •• Average collective control (5) 
Difference Variables: 

Ax n : Horizontal separation (6) 

Az « Zfl, - Zg = Vertical separation . (7) 

AS = e m - 3 S : Differential pitch attitude . (8) 
A8 lc a fltem" Pics : Differential cyclic control (9) 
A0 C * 6^- 0^: Differential collective control (10) 
where |J5 j c |< IS degrees and jG c |< 10 degrees 
for each helicopter. 
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Fig. 2: Model for TUtS longitudinal dj-nanics. 

The variables defined above are used to derive the 
equations which govern the planar motion of the 
TIHS. In doing so it is assumed that tne trim 
condition is one in whjch the helicopters are near 
hover with the tethers vertical and the spreader bar 
horizontal. It is also assumed that the helicopter 
cg.*s lia on the main rotor drivushafts. *his 
assumption forces the helicopters to h3ve rw nominal 
pitch for the trim con ition described above. With 
these assumptions a linear state space model is ' 
constructed having the following form; j 

^p = * Bp lip : SpfR 12 . Up£R 4 01) . 

c[2z SAx A0 Ax A9 BEO AZ E£ Ex E9 Az E£l T (12) 
Up a l EO c a A8 lc a A0 C E3j c 1 T (13) 

Ap« diagfApj.Ap^Ap,;) (H) ; 

Bp “ ciag (Op j , Dp 2 , Bpj). (IS) 

Note that due to the symmetry or the "equal tether 
c onfiguration* . the modal decouples into three basic 
subsystems to be discussed subsequently. It should 
also be noted that the components of j$p represent 

small perturbations from nominal trim values; i.e. 
the values when the tethers arc vertical and the 
spreader bar horizontal. 





72 Three Basic Subsustems and Output S election 
The three basic subsystems of the TLHS are now 
discussed, and variables to be commanded by the 
pilot (outputs) are selected. 

The pair (Ap ( . Bp]) describes the system's •average 

vertical motion* (AVtl). This mcticn involves the Ez 
degree of freedom, and is controlled by issuing 
average collective commands (E9 C ). The subsystem 

thus has ona input which can be used to control at 
most one output, Ez. This SISO first order subsystem 
is stable and characterizes the natural damping that 
occurs during vertical climbs. 

The pair (Ap^ B^) describes the. system’s 
"summetric motion* (Stt) [lL This motion involves tho 
Ax and A9 degrees of freedom, and is controlled by 
issuing differential cyclic commands (Aglc). The 
subsystem thus has only ona input which can ba used 
to control at most one output. Tha variable Ax is 
selected" as tha output because maintaining 
horizontal separation is critical. This - SISO fourth 
order subsystem is unstable and is characteristic of 
any helicopter tethered to a fixed point in space. The 
instability is referred to as the “ tethered helicopte r 
mode* (Fig. 3). This mode is a result of the tethers 
being attached a distance h below the helicopter 
cg.'s (Fig. 2) [2L 
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Fig. 5: Syctetric Motion Foies S Zeros. 

Finally, the pair (Ap 3 . B p3 ) describes the system's 
*3nti-summgtric motio n* (ASn). This motion irn’olves 
the Ex. £9. Az. Zi degrees of freedom, and is 
controlled by issuing differential collective 
commands (A9 C ) and average cyclic commands 

(EP tc ). The subsystem thus has two inputs which can 

be used to control at most two outputs. The variables 
x^-Ex and Ex are selected as the outputs because 
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they are fundamentally related to the problem of 
controlling master/slave workload during horizontal 
flight. The quantity x L *Ix represents- the Ussl 

deviation trom the center and is given by: 
x L - Ex = h E9 * (Z/L) Az»HE£ (IG) 

where h. 2. L and H are constant geometric 
parameters (Fig. 2). This two-input two-output 
(TITO) seventh order subsystem is unstable, and if 
Z = 0, is characteristic of any single hovering 
helicopter carrying a sling load. The instability hero 
is referred to as the * backfieppinn mode* (Fig. 4) 
since it is due to the backflapping of the helicopter 
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The outputs described above are written in matrix 
ferm as follows: 

Up B CpXp; Up£R 4 . _ (17) 

where ^ » t Ez a Ax a x|_ - Ex Ex V . 

Cp * dbg (Cpj, Cp 2 . Cp 3 ) (19) 

and Cp,. Cp 2 , Cp 3 are associated with the AVtl, Stt, 

and ASn subsystems, respectively. 

Given the above state space model (Ap. Bp, Cp), the 

TLHS (plant) transfer fur. .ion matrix is then defined 
as follows-" 

G p ■ Cp (si - Ap ) 1 Bp (20) 

= diag (G ? , (s). G p2 (s). G p3 (s)) (2 1 ) 

where the (G p j(s), i= 1,2,3) specify the input/cutput 

characteristics of the three subsystems discussed 
above, and are given by: 


G pi E ^pi^ 5 * " 


I* 1.2. 3. 


It is emphasized at this point that the above 
decoupling into two 5150 systems and one TITO 
system occurs because the helicopters are assumed 
to be identical and the tethers arc assumed to have 
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the same length. The ' unequal tether problem* . in 
whicn the Stl and ASM couple to form a three-input 
three-output subsystem. is addressed in [41. 

It is strongly emphasized that although the equations 
are decoupled in terms of the defined variables, they 
remain coupled in certain quantities such as the 
individual helicopter collective and cyclic controls. 
Consequently all three subsystems and their 
corresponding compensators must be evaluated 
simultaneously so that important quantities can be 
•untangled* using eqs. (1) - (10). This coupling, in 
effect, makes controlling the longitudinal twin lift 
’equal tether configuration' a four-inout fcur-output 
design problem. 

3. AFCS Structure and Design Specifications 
3.1 AFCS Structure and Function Definitions 
The final control system structure is shown in Fig. 5. 
In this liino negative feedback structure there is the 
plant Gp(s) (TLHS), a dynamic augmentation 

consisting of four integrators (one per channel), and a 
dynamic 16 th order compensator Klqg/LTR^* 



The compensator is designed sc that the closed loop 
system meeis the specifications- presented below. 
Since there are three subsystems, the LQG/LTR 
compensator has the fc.m: 

*'LQG/LTr( s ^ = Ofj(s)» ^(^ (23) 

where the Kj(s) denote the individual subsystem 
compensators. 

The TLHS (plant) plus the integrators is referred to 
as the design plant, its transfer function matrix is 


defined as follows: 

G(s) a C(s! - A)“'B (24) 

= diag (Gj(s), G 2 (s), G 3 (s)) (25) 

* Gp(s)/s (26) 

where A » diag (Aj. A 2 . A 3 ) (27) 

B = diag (Bj. B 2 . B 3 ) (2B) 

C ■ diag (C,. C 2 . C 3 ) (29) 


and the design plan t input/output characteristic 
is given by- 


Gj(s) a Cj(sl - Aj)”^Bj (30) 

= Gpj(s)/s i =1. 2. 3 (31) 

where i - l. 2. 3 denote the AVM, sn. and ASM 
subsystems respectively. 

In an analogous manner, the final loop, sensitivity, 
and closed loop transfer function matrices can ba 
defined as follows: 

Gj_(s) a diag (G^j(s). Gj^s). G^ 3 (s)) (32) 

Sp(s) « diag (S F , (s). S^s). S^fs)) (33) 

Tp(s) * diag (Tjq(s), Tj^fe). Tpnj(s)). (34) 

32 Performance Specifications - 
To develop reasonable performance specifications. 
i.e. specifications which are 'passenger friendly', the 
plant singular values (Fig. 6) were analyzed and 
closed loop transient responses were studied to get a 
feel for typical twin lift pilot maneuvers. 
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0 ns performance specification is zero steady state 
error to pilot step reference commands, r, for £z. Ax, 
x^- lx. and lx. Another is complete rejection of step 

disturbances, <L with respect to the above four output 
variables. The four integrators introduced in the 
AFCS (one per charnel) guarantee that both of these 
performance specifications are met (see Internal 
Model Principle [7D. The integrators also help with 
low frequency command following, low frequency 
disturbance rejection, and low frequency sensitivity. 
To guarantee that these properties are built into the 
design, it is required that the final subsystem 
sensitivity functions satisfy the following frequency 
domain constraints: 

|5pj(ju)! < -20 db for all u> < .05 rad/sec (35) 

| s F2 (ja>) J < -20 db for all o < .04 rad/sec (36) 
C max l 5 (: 3 Cj{ii)l < -20 db for all u> < .06 rad/sec (37) 



where cJ max [ • 1 denotes the maximum singular value. 

3.2 Rcoustness Specifications 

To tie robust with respect to the high frequency 

unmodeled main rotor dynamics (at o> ~ 7.7 rad/sec), 

the final subsystem closed loops must satisfy the 

fallowing frequency domain constraints: 

jTpjCjtia) | < -20 db for all ta > 10 rad/sec (38) 

|Tp 2 (jtjj)| < -40 db for all u > 10 rad/sec (39) 

cJmaxf r F 3 (i^)l < *20 db for all tu > 10 rad/sec (40) 
These also help in attenuating high frequency sensor 
noise, a 

The above ptrformance and stability robustness 
specifications are nonsensical unless the nominal 
closed loop system is stable . This . however, is not 
a problem since the LQG/ITR design methodology 
will be used. The methodology guarantees nominal 
stability and offers excellent stability margins (7J. 

4. LQG/LTR Design Methodology 
4.5 Obtainioo the design 

Next the LQG/UTR design methodology is applied to 
the minimum phase [61 [91 [l 01 design plant G(s) = 
C(s l- a)” ' 3 (i.e., TIHS plus integrators}:- The design 
process consists of two steps: 

Stop 1- Design a 'nice target loop', G^s) = C(sl-A)H, 

where H is found by solving a filter algebraic 
Riccati equation. The Riccati equation depends on 
design parameters that can be manipulated to place 
loop transmission zeros at appropriate locations to 
give us loop shapes that satisfy the aforementioned 
'passenger friendly* specifications. Gj^p(s) has 

guaranteed closed loop robustness properties. 

Steo 2: "Recover* the nice loop properties of Gj^s) 

from the design plant: 

LIB- G(s)K lqg/ltr (s) > G kf (s) (41) 

as 

Pi* P2* ?3 >0: 

i.e. as the design parameters P 1 .P 2 .P 3 (on® fcr 
each subsystem) get sma ll th e loop transfer function 
matrix approaches a function with good performance 
and guaranteed robustness properties. 

4.2 Evaluation of the final design 
The final loop, sensitivity, and closed loop singular 
values appear in Figs. 7a - 7c The shaded regions 
indicate that the design specifications (eqs. 
(37)-(40)) have been met. 

Fig. 7a shows that all loop singular values have 
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gain-crcsscver frequencies less than t.S rad/sec and 
above .6 rad/sec. Their low frequency slopes of -20 
db/dec confirms the presence of the four integrators 
in the loop. . 
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Fig. 7a: Final loop Sing u lar Va lues. 

Fig. 7b shows that all sensitivity singular values lie 
below the -20 db ' performance region* thus 
guaranteeing good low frequency command following, 
good low frequency disturbance rejection, and good 
low frequency sensitivity reduction. The pox 511 and 
ASM sensitivity properties at frequencies above .6 
rad/sec are expected since these subsystems have 
unstable modes which force the tradeoff cf nice 
feedback _ properties at high frequencies fcr nice 
feedback properties at low frequencies [91 [IOL 



_ equeney (rad/sec) 

Fig. 7b: Fir.il Sensitivity Singular Values; 

Fig. 7c shows that all closed loop singular values lie 
below the -20 db 'robustness reoion* thus 
guaranteeing robustness to the unmodeled rotx 
dynamics st u ~ 27 rad/sea The design also provides 
a degree of high frequency sensor noise attenuation 
built into the design. 

To complete the evaluation of the design, the 
following reference commands were issued, assuming 
zero initial conditions, no disturbances (cl = 0 ). and 
no sensor noise (q = Q): 
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_Flg, 7c: Final Closed Loop Singular Values. 
Reference Command 
ASflt Ez = 5 ft/sec 

S0. : Ax = 2 ft 


ASH: XL-Ex = -lft 

Ex = 5 fl/sec 

The results of the above linear simulation are given 
in Figs. 8a - Oc and 9a - 9d. It is emphasized that 
cqs. (l) - (10) are used to ’untangle* the variables of 
interest (e.g. the individual helicopter velocities and 
controls). 

initially the helicopters are at rest with the tethers 
vertical and the spreader bar horizontal. The pilot 
reference commands being issued ask for an effective 
system steady state motion of 7.071 ft/sec in a 
direction 45 degrees with respect to the horizontal, 
with an incremental horizontal separation cf 2 feet 
(Ax = 2). and the payload shifted 1 foot to the left of 

the center^ - Ex = -l). _ ... 

Figures 3a and 9a show how the AFCS coordinates the 
helicopter controls so as to achieve the outputs 
commanded by the pilot The initial control 
transients, although tolerable, are relatively large but 
this is expected since the AFCS is trying to achieve a 
vertical climb of 5 ft/sec, a horizontal drive of 5 
ft/sec. keep the load 1 ft behind the center, all while 
maintaining an incremental horizontal separation of 2 
feet between two inherently unstable 14000 lb 


Dlackhawksl 

Fig. Qa shows that the AFCS keeps the master and 
slave collectives out of phase throughout the 10 
second maneuver. Fig. 8b shows the same for the 
helicopter velocities. 

The ret effect is seen in Fig. 8c which shows how the 
AFCS effectively uses the Az degree of freedom to 
control the Ex load motion. After approximately 6 
seconds, the Az and x^- Lx motions come into phase 




rig. 8b: Master 5 Slave Vertical Velocities. 



with one another and decay to their correct steady 
state values of -2 feet anJ -1 foot, respectively. The 
motion observed is essentially the vertical spring 
mode associated with the anti-symmetric motion 
(Fig. 4). In the steady state a difference in collective 
controls provides the force to maintain the slave 
above the master (Az = -2). This in turn maintains the 
payload to the left of the center (x^Ix = -1). The key 

point to be made here is that the AFCS must 
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•tradeoff* tne system £2 motion in order to follow 
pilot reference commands in ^ - Ex. It truis follows 

intuitively that relatively poor vertical acceleration 
characteristics may result as a consequence of 
requiring the use of vertical separation to control the 
load deviation from the center. This is somewhat 
apparent in figure 8b which exhibits the expected 
vertical acceleration characteristics. It must be 
pointed out that the vertical acceleration 
characteristics can be improved by increasing the 
bandwidth of the rinal AVH and ASH loop, but only at 
the expense or larger collective controls, and the 
possible excitation of the high frequency unmodeled 
rotor dynamics. 

Fig. 9a shows that the AFCS initially maintains the 
helicopter cyclics slightly out of phase. In doing so 
both helicopters pitch downward (Fig. 9b) and 
accelerate forward (Fig. 9c). In trying to achieve a 2 
foot horizontal separation, the AFCS coordinates the 
cyclics so as to induce the necessary helicopoter 
pitching motions (Fig. 9c). Due to the inertia of the 
helicopters, however, a “horizontal separation 
overshoot* occurs (Fig. 9d). This overshoot can be 
decreased substantially by increasing the bandwidth 
of the final Sfl but this can be done only at the 
expense of substantially larger cyclic controls and the 
possible excitation of the high frequency unmodclsd 
rotor dy namics. 
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Fig. 9a: Master 5 Slave Cyclic Controls. 

Figs. 9a - 9d show that the AFCS must essentially 
•tradeoff the pitching motions of the helicopters in 
order to follow pilot reference commands for Ax. 
8ecause of this tradeoff, it follows that a sluggish 
horizontal acceleration characterictic may result for 
a low bandwidth SM design (Fig. 9c). As in the vertical 
case discussed above, it should be noted that the 
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Fig. 9c: Master 5 Slave Horizontal Velocities. 
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Fig.Sd: Horizontal Separation. 

horizontal "acceleration charateristics may be greatly 
improved by increasing the bandwidth of the Sfl loop, 
but only at the expense of larger cyclic controls, 
larger pitch rates, and the possible excitation of the 
high frequency unmodeled main rotor dynamics. 
Comparing figures 8b & 9c, it is seen that the vertical 
acceleration charatcristics are slightly worse than 
the horizontal acceleration charateristics. This is 
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primarily due to the fact that although Az does appear 
in the expression for Zx (eq. 16). it is not 

directly ’penalized*; i.e. it is not an output as is its 
horizontal counterpart. Ax. 

It should he noted that the four integrators (one per 
input channel) have provided the desired zero steady 
state error to the pilot's step reference commands. It 
should also be noted that because of the unstable 
’tethered helicopter* and "backflapping* modes, a 
minimum bandwidth is required to just stabilize the 
’symmetric and anti-symmetric motions. This, of 
course, translates into minimum control rates and 
more importantly (from a passenger's viewpoint), 
larger pitch rates. These large pitch rates must be 
lived with, and are linear functions of the pilot 
reference commands as well as the system initial 
conditions. 

To improve the vertical performance one might try 
first to use a configuration in which the slave tether 
is longer than the master tether (41. Such a 
configuration greatly improves the relative stability 
of the nominal model by introducing damping in the 
’tethered helicopter* mode and the “backflapping* 
mode, as well as the other oscillatory modes of the 
TLHS. It must be emphasized that by bringing' the 
unstable poles of the TLHS closer to the left half 
plane, the minimum bandwidth needed to stabilize our 
system decreases. This is desirable because a lower 
bandwidth implies smaller controls and hence smaller 
pitch- rates. 

Another great improvement would result by flying, 
with the helicopter-tether attachment points directly 
at the Blackhawk cg.’s (h = 0). This has been proposed 
at Sikorsky Aircraft and is dicussed in 111. . 

Another idea for improving the performance of the 
TLHS is to use real lime tether length control. This 
would introduce another degree of freedom in the 
controls which v/ould seem to be of tremendous help 
in controlling the payload-bar assembly. 

S. Summary 

This paper has shown that the LQG/LTR design 
methodology, coupled with singular value ideas, can 
be used to formulate performance and stability 
robustness specifications which can then be used to 
develop a centralized nino AFCS for a twin lift 
hchicopter system (TLHS). The methodology is 
extremely powerful in that it possesses enough 
degrees of freedom to tailor designs to prescibed 


pilot specifications. The paper and (4l show that the 

methodology can be used to address the TLHS control 
problem with repect to the following five 
fundamental feedback issues: 

(1) low frequency command following. 

(2) low frequency disturbance rejection, 

(3) low frequency sensitivity reduction, 

(4) robustness to high frequency modeling, and 

(5) attenuation of hign frequency sensor noise. 
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